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ABSTRACT 

We present precise Doppler measurements of four stars obtained during the past decade at Keck 
Observatory by the Cahfornia Planet Survey (CPS). These stars, namely, HD 34445, HD 126614, 
HD 13931, and G1179, all show evidence for a single planet in Keplerian motion. We also present 
Doppler measurements from the Hobby-Eberly Telescope (HET) for two of the stars, HD 34445 and 
Gl 179, that confirm the Keck detections and significantly refine the orbital parameters. These planets 
add to the statistical properties of giant planets orbiting near or beyond the ice line, and merit follow- 
up by astrometry, imaging, and space-borne spectroscopy. Their orbital parameters span wide ranges 
of planetary minimum mass (M sini = 0.38-1. 9Mjup), orbital period (P = 2.87-11.5yr), semi-major 
axis (a = 2.1-5.2 AU), and eccentricity (e = 0.02-0.41). HD 34445b (P = 2.87yr, M sin i = 0.79 Mj^p, 
e = 0.27) is a massive planet orbiting an old, G-type star. We announce a planet, HD 126614 Ab, and an 
M dwarf, HD 126614 B, orbiting the metal-rich star HD 126614 (which we now refer to as HD 126614 A). 
The planet, HD 126614 Ab, has minimum mass Msini — 0.38 Mjup and orbits the stellar primary with 
period P = 3.41yr and orbital separation a = 2.3AU. The faint M dwarf companion, HD 126614 B, is 
separated from the stellar primary by 489 mas (33 AU) and was discovered with direct observations 
using adaptive optics and the PHARO camera at Palomar Observatory. The stellar primary in this 
new system, HD 126614 A, has the highest measured metallicity ([Fe/H] = -1-0.56) of any known planet- 
bearing star. HD 13931 b (P = 11.5 yr, M ami = 1.88 Mjup, e = 0.02) is a Jupiter analog orbiting a near 
solar twin. G1179b (P = 6.3yr, Af sin i = 0.82 Mjup, e = 0.21) is a massive planet orbiting a faint M 
dwarf. The high metallicity of Gl 179 is consistent with the planet- metallicity correlation among M 
dwarfs, as documented recently by Johnson & Apps. 

Subject headings: planetary systems ~ stars: individual (HD 34445, HD 126614, HD 13931, G1179) 
— binaries: visual — techniques: radial velocity — techniques: photometric — 
techniques: high angular resolution 



1. INTRODUCTION 

The distributions of the masses and orbits of jovian- 
mass exoplanets offer key tests of planet formation the- 
ory. Most theories predict that giant planets form be- 
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yond the "snow line" and migrate inward on a time scale 
that competes with the lifetime of the protoplanetary 
disk (Thommes et al. 2008; Ida & Lin 2008; Rice & Ar- 
mitage 2005; Alibert et al. 2005; Trilfing et al. 2002). 
Among various theories for formation, core accretion has 
been shown efficient at producing, within ~3 Myr, plan- 
ets of Neptune to Jupiter mass, orbiting within 5 AU 
(Benz et al. 2008; Dodson-Robinson et al. 2008). 

These models of giant planet formation and orbital 
evolution may be directly tested against observations 
of giant planets found by the Doppler method. The 
models predict a clearing of gaps in the protoplanetary 
disks, establishing the distribution of planet masses, al- 
lowing a direct comparison with Doppler observed min- 
imum masses (Afsini). According to the models, the 
disk dissipates after the planets have undergone some 
inward migration, leaving them at their current orbital 
distances. The resulting predicted distribution of semi- 
major axes can be compared with the observed orbits 
of giant planets. There remains potential value in both 
enhancing the sophistication of planet formation theory 
and in observing a large enough statistical sample of gi- 
ant planets to permit robust and informative tests of the 
theory. Moreover, planet-planet interactions among gi- 
ant planets must be predicted and compared with the 
distributions of orbital elements (especially eccentricity) 
for systems containing multiple giant planets, e.g. Wright 
et al. (2009); Ford & Chiang (2007); Ford & Rasio (2008); 
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Juric & Treniaine (2008). Giant planets also gravitation- 
ally interact with the dust in their planetary system to 
shape, on timescales of only years, the dust evolution of 
the planetary system (Lisse et al. 2007; Beichman et al. 
2007; Payne et al. 2009). 

As of May 2009, 350 exoplanets have been discov- 
ered, with remarkable properties including close-in or- 
bits, large orbital eccentricities, multi-planet systems, 
and orbital resonances (Mayor & Udry 2008; Marcy et al. 
2008). The hot jupiters have received the most attention, 
observationally and theoretically, yielding extraordinary 
information about their chemical composition, internal 
structure, atmospheric behavior. However most known 
gas giant planets orbit beyond 1 AU, realizing the pop- 
ulation that formed beyond the ice line as predicted by 
theory, and as seen in our Solar System. 

In 1997, we began a Doppler search for giant planets 
in Jovian orbits at Keck Observatory. We monitor over 
1800 stars within 50 pc, with special attention given to 
those within 20 pc. We have acquired a Doppler time 
baseline of well over Syr for nearly all of them. The 
detected long-period exoplanets reveal the distribution 
of their masses, semimajor axes, and orbital eccentrici- 
ties for the general population of planetary systems. Re- 
markably, the exoplanets exhibit a sharp rise in occur- 
rence beyond 1 AU (Johnson et al. 2007b; Gumming et al. 
2008), indicating a great population of giant planets re- 
sides there. Many of these planets remain undiscovered 
even after 10 yr of Doppler monitoring, as the amplitudes 
of a few meters per second require high Doppler preci- 
sion and a clear indication of Keplerian motion, which is 
challenging for orbital periods comparable to the dura- 
tion of observations. Nonetheless, the analysis of Doppler 
completeness shows that 15-18% of all nearby stars have 
giant planets between 3-20 AU (Gumming et al. 2008). 
Surely, these giant planets offer strong statistical infor- 
mation on the formation and subsequent dynamical evo- 
lution of gas giants in general. 

Unfortunately, with only a decade of data orbits be- 
yond 4AU are just coming into our Doppler field of 
view. The recently announced Jupiter-analog orbiting 
HD 154345 with a = 5.0AU (and a circular orbit) exhib- 
ited nearly one full orbital period only after 10 full years 
of Doppler data were collected (Wright et al. 2008). But 
the number of giant planets known orbiting beyond 1 AU 
remains so small that extraordinary statistical efforts are 
required to extrapolate the true underlying properties 
(Gumming et al. 2008). Thus, there remains a need 
for enlarging the observed population of giant planets, 
especially with the prospect of follow-up observations 
by such instruments as the James Webb Space Tele- 
scope (JWST), the Gemini Planet Imager (GPI), and 
the Spcctro-Polarimctric High-contrast Exoplanct Re- 
search (SPHERE) instrument for the Very Large Tele- 
scope (VLT). 

In the future, knowledge of giant planets around the 
nearest stars will be crucial for detecting Earth-sized 
planets, e.g. by the Space Interferometry Mission (Un- 
win et al. 2008), as the giant planets add a "noise" to 
the astrometric signal (Traub et al. 2009). One concern 
is that multiple giant planets orbiting beyond 2 AU will 
cause curved astrometric motion, with the linear part be- 
ing (erroneously) absorbed into the inferred proper mo- 
tion. The resulting residuals will have a time scale of 



^lyr, constituting an astrophysical "noise", compromis- 
ing the detection of the terrestrial planets. Thus, future 
astrometric missions will benefit greatly from ~15yr of 
radial velocity (RV) data. The characterization of both 
giant and rocky planets will be important for future mis- 
sions that image and take spectra of planets, such as 
the Terrestrial Planet Finder and Darwin (Kaltenegger 
et al. 2006; Lawson et al. 2008). As a result we con- 
tinue to survey nearby stars that are likely targets for 
such surveys (Kaltenegger et al. 2010). Here we describe 
Doppler measurements from Keck Observatory and from 
the Hobby-Eberly Telescope (HET) for four stars that 
shows signs of harboring a planet beyond 1 AU. 

2. STELLAR SAMPLE AND PROPERTIES 

Among the 1800 FGKM stars we monitor at the Keck 
telescope for giant planets, 1330 of them have been mon- 
itored from 1997 to the present with a precision of 1- 

3 ms^^. The sample contains a nearly volume-limited 
sample of stars within five bins oi B — V color, between 
0.60 and 1.5. It is nearly devoid of both binary stars sep- 
arated by less than 2" due to mutual contamination at 
the slit of the spectrometer and magnetically active stars 
younger than 1 Gyr due to their excessive velocity "jit- 
ter" (Wright 2005). The complete list of target stars is 
given in Wright et al. (2004) including the stars that do 
not have detected planets, permitting statistical analyses 
of the frequency of planets, e.g. Johnson et al. (2007a); 
Gumming et al. (2008); Butler et al. (2006). 

Only 7% of our stars surveyed for a decade reveal clear 
Keplerian Doppler variations at a 3-a limit of 10 ms^^. 
However Gumming et al. (2008) accounted for incom- 
pleteness in our Doppler survey and predicted that -^18% 
of FGK stars harbor a gas giant of saturn mass or above 
within 20 AU. This prediction suggests that another 11% 
of our target stars harbor giant planets yet to be revealed, 
with continued monitoring. The four stars presented here 
presumably harbor some of these giant planets by pre- 
dicted Gumming et al. (2008). 

We measure atmospheric parameters of the target stars 
by LTE spectroscopic analysis of our Keck/HIRES spec- 
tra using the "SME" (Spectroscopy Made Easy) code 
(Valenti & Piskunov 1996) as implemented by Valenti & 
Fischer (2005) and Fischer & Valenti (2005). 'The analy- 
sis yields a best-fit estimate of Tcff, logg, [Fe/H], vsini. 
We could not carry out an LTE analysis of Gl 179 because 
its cool temperature resides outside the lowest temper- 
ature for which our continuous and molecular opacities 
are accurate. 

The luminosity of each star is determined from the 
apparent V-band magnitude, the bolometric correction, 
and the parallax from Hipparcos. Prom Tge and the lumi- 
nosity, we determine the stellar mass, radius, and an age 
estimate by associating those observed properties with 
a model from the stellar interior calculations of Takeda 
et al. (2007, 2008). These properties, along with paral- 
laxes and implied distances from the van Leeuwen (2007) 
reduction of Hipparcos data, are tabulated in Table 1. 

Since G1179 is faint and cool, we could not use 
SME or the Takeda et al. stellar interior models to 
compute stellar properties. Instead, we used a va- 
riety of other techniques to estimate the properties 
listed in Table 1. We determined the mass of G1179 
by using the mass-luminosity calibration of Delfosse 
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TABLE 1 
Stellar Properties 



Parameter 


HD 34445 


HD 126614 A 


HD 13931 


Gl 179 




GO 


KO 


GO 


M3.5 


Mv 


4.04 ± 0.10 


4.64 ± 0.17 


4.32 ± 0.10 


11.5 ± 0.11 


B — V 


0.661 ±0.015 


0.810 ± 0.004 


0.642 ± 0.015 


1.590 ± 0.015 


V 


7.31 ± 0.03 


8.81 ± 0.002 


7.61 ± 0.03 


11.96 ± 0.03 


Parallax (mas) 


21.5 ±0.7 


13.8 ± 1.0 


22.6±0.7 


81.4±4.0 


Distance (pc) 


46.5± 1.5 


72.4 ±5.3 


44.2 ±1.4 


12.3±0.6 


[Fc/H] 


±0.14 ±0.04 


+0.56 ±0.04 


+0.03 ± 0.04 


+0..30±0.10 


Teff (K) 


5836 ±44 


5585 ± 44 


5829 ± 44 


3370 ± 100 


V sin i (kms 


2.7±0.5 


2.0 ±0.5 


2.0 ±0.5 


<1.5 


logg 


4.21 ±0.08 


4.39 ±0.08 


4.30 ±0.08 


4.83 


L. (Lq) 


2.01 ±0.2 


1.21 ±0.19 


1.57 ±0.14 


0.016 ±0.02 


(Mq) 


1.07 ±0.02 


1.145 ±0.03 


1.022 ±0.02 


0.357 ±0.03 


R* (Rq) 


1.38 ±0.08 


1.09 ±0.06 


1.23 ±0.06 


0.38 ±0.02 


Age (Gyr) 


8.5 ±2.0 


7.2 ±2.0 


8.4 ±2.0 




logi?' 
Prot (days) 


0.148 


0.152 


0.161 


0.956 


-5.07 


-5.44 


-4.99 


-5.20 


~22 


~99 


~26 





et al. (2000), applying its apparent K-band magnitude, 
K = 6.942, and parallax of 81.4 mas. The resulting mass 
is Mi, = 0.357 ± 0.030 Mq, with the uncertainty from the 
scatter in best-fit M-L relation. Johnson & Apps (2009) 
estimate a metallicity of [Fe/H] = 0.30 ± 0.10 for G1179 
based on its absolute K-band magnitude, Mk, and 
V-K color. We estimate = 0.016 ± 0.002 Lq from the 
bolometric correction, Teg = 3370 ± 100 K from Bessell 
(1995), i?* = 0.38 ± 0.02 i?o from the Stcfan-Boltzmann 
law, iisini < 1.5kms~^ from visual inspection of the 
G1179 spectrum (Figure 1), and log fif = 4.83 from the 
stellar mass and radius. 

We also measure the chromospheric emission in the 
Ca II H & K line cores, providing S'hk values (Isaacson 
2009; Wright et al. 2004) on the Mt. Wilson system. For 
each star, the time series Keck RVs and Shk values are 
uncorrelated (as measured by a Pearson correlation co- 
efficient). We converted the average 5hk values for each 
star to logi?^K per Noyes et al. (1984), providing an 
estimate of the age and rotation period of the stars. The 
log relation is not calibrated for cool M stars, but we 
include the computed value for Gl 179 nonetheless. All of 
the resulting stellar parameters are reported in Table 1. 
We discuss the salient properties of each star in § 5-8. 

HD 126614 is is a binary star system composed of 
a bright primary, HD 126614 A. and a faint M dwarf 
companion, HD 126614 B, separated by ^0'.'5. The 
heretofore unknown companion is significantly fainter 
{AV = 7.8 mag) and did not significantly contaminate 
the spectroscopic or photometric observations. (See §6 
for details.) 

3. PHOTOMETRIC OBSERVATIONS FROM FAIRBORN 
OBSERVATORY 

We have obtained between two and seven years of high- 
precision differential photometry of all target stars in this 
paper with the exception of Gl 179. The observations 
were acquired with the T12 0.8 m automatic photomet- 
ric telescope (APT) at Fairborn Observatory. This APT 
can detect short-term, low-amplittidc brightness variabil- 
ity in solar-type stars due to rotational modulation in 
the visibility of photospheric starspots (e.g., Henry et al. 
1995) as well as longer-term variations associated with 
stellar magnetic cycles (Henry 1999). Therefore, pho- 
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Fig. 1. — Spectra near the Ca II H line for all four stars discussed 
here. The emission reversals reflect magnetic activity on stars that 
correlates with photospheric velocity jitter. The resulting chromo- 
spheric cooling rate, measured as a fraction of stellar luminosity, 
log Rjjj^, is listed in Table 1. The inferred values of jitter are listed 
in Table 3. 

tometric observations can help to establish whether ob- 
served radial velocity variations in a star are due to reflex 
motion caused by a planetary companion or due to the 
effects of stellar activity (e.g., Queloz et al. 2001; Paulson 
et al. 2004). Photometric observations can also lead to 
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TABLE 2 

Summary of Photometric Observations from Fairborn Observatory 



Target 


Comparison 


Comparison 






Date Range 


Target cTghort 


Comp CTshort 


Target CTiong 


Comp CTiong 


Star 


Star 1 


Star 2 


A^obs 


A^yrs 


(JD-2,400,000) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


HD 34445 


HD 34907 


HD 35241 


560 


7 


52591-54913 


0.00155 


0.00156 


0.00102 


0.00120 


HD 126614 


HD 127265 


HD 124988 


113 


2 


54171-54634 


0.00158 


0.00163 


0.00002 


0.00161 


HD 13931 


HD 13024 


HD 14064 


181 


2 


54730-54532 


0.00153 


0.00164 


0.00158 


0.00060 



the detection of planetary transits and the direct deter- 
mination of planetary radii, as in Winn et al. (2008). 

We acquired 854 nightly photometric observations of 
the three planetary host stars between 2000 November 
and 2009 March. HD 13931, HD 34445, and HD 126614 A 
were observed by the T12 APT, a functional duplicate 
of the T8 0.8 m APT and two-channel photometer de- 
scribed in Henry (1999). All target stars were observed 
differentially with respect to two comparison stars in the 
Stromgren b and y bands. One to four observations per 
night were made of each target, depending on the sea- 
sonal visibility of each star. To increase our photometric 
precision, we averaged the individual b and y differential 
magnitudes into a composite {b + y)/2 passband. 

The results of our APT observations are summarized in 
Table 2, which identifies the comparison stars used, the 
number and timespan of the observations, and our mea- 
sures of brightness variability in the target and compari- 
son stars. AH three of the observed stars have very little 
chromospheric activity (Table 1 and Figure 1). There- 
fore, given the well-known connection between chromo- 
spheric activity and brightness variability in solar-type 
stars (see e.g., Henry 1999, Fig. 11), we did not expect 
to find significant photometric variability in these four 
targets. 

We use the standard deviation as a simple metric for 
brightness variability of our target and comparison stars, 
and we measure the variability on both night-to-night 
(cshort) and year-to-year (criong) timescales. In general, 
the standard deviations of individual measurements of 
pairs of constant stars fall in the range 0.0012-0.0017 
mag for these telescopes. Column 7 of Table 2 gives 
the standard deviation of the individual target minus 
comparison star differential magnitudes averaged across 
both comparison stars and across all observing seasons. 
Column 8 gives the standard deviation of the individ- 
ual differential magnitudes between the two comparison 
stars averaged across all observing seasons. We com- 
pute the (cTiong) values in the same manner but with the 
yearly mean differential magnitudes rather than individ- 
ual nightly magnitudes. The photometric results for each 
of the four stars are briefly discussed in the following sec- 
tions. Given the moderately long orbital periods and the 
current uncertainties in the orbital parameters (Table 3), 
a search for planetary transits in any of the four target 
stars is impractical. 

4. SPECTROSCOPIC OBSERVATIONS AND KEPLERIAN 
ORBITAL SOLUTIONS 

4.1. Keck Observations and Doppler Reduction 

We observed HD 34445, HD 126614 A, HD 13931, and 
Gl 179 using the HIRES echelle spectrometer (Vogt et al. 
1994) on the 10-m Keck I telescope. These stars were 



each observed 30-70 times each over 10-12 yr. All obser- 
vations were made with an iodine cell mounted directly 
in front of the spectrometer entrance slit. The dense 
set of molecular absorption lines imprinted on the stel- 
lar spectra provide a robust wavelength fiducial against 
which Doppler shifts are measured, as well as strong con- 
straints on the shape of the spectrometer instrumental 
profile at the time of each observation (Marcy & Butler 
1992; Valenti et al. 1995). 

We measured the Doppler shift from each star-times- 
iodine spectrum using a modeling procedure modified 
from the method described by Butler et al. (1996). The 
most significant modification is the way we model the in- 
trinsic stellar spectrum, which serves as a reference point 
for the relative Doppler shift measurements for each ob- 
servation. Butler et al. use a version of the Jansson 
(1995) deconvolution algorithm to remove the spectrom- 
eter's instrumental profile from an iodine-free template 
spectrum. We instead use a new deconvolution algo- 
rithm developed by one of us (J.A.J.) that employs a 
more effective regularization scheme, which results in sig- 
nificantly less noise amplification and improved Doppler 
precision. 

Figure 2 shows the radial velocity time series for four 
stable stars with characteristics similar to the planet- 
bearing stars presented herein, demonstrating our mea- 
surement precision over the past decade on a variety of 
spectral types. In 2004 August, the Keck HIRES spec- 
trometer was upgraded with a new detector. The pre- 
vious 2K X 2K pixel Tektronix CCD was replaced by 
an array of three 4K x 2K pixel MIT-LL CCDs. The 
new detector produces significantly higher velocity pre- 
cision due to its improved charge transfer efficiency and 
charge diffusion characteristics, smaller pixels (15 fim 
vs. 24 /im), higher quantum efficiency, increased spec- 
tral coverage, and lower read noise. Our post-upgrade 
measurements exhibit a typical long-term rms scatter 
of ~1.5ms^^ for bright, quiescent stars, compared to 
~2.5ms~^ for pre-upgrade measurements. (Measure- 
ments prior to JD 2,453,237 are pre-upgrade.) The pre- 
and post-upgrade measurements also lack a common ve- 
locity zero point, but we fit for and corrected this offset 
to within ~2ms^^ for every star observed at Keck us- 
ing a large set of stable stars with many pre- and post- 
upgrade observations. To further limit the impact of the 
velocity discontinuity, we let the offset float in the Ke- 
plerian fits below, effectively treating pre-upgrade and 
post-upgrade observations as coming from two different 
telescopes (except for HD 13931, which does not have 
enough post-upgrade measurements to allow for a float- 
ing offset). 

With the exception of G1179, these planet-bearing 
stars are bright G and K stars, and the Keck RVs 
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Fig. 2. — Radial velocity time series for four stable stars in our 
Keck Doppler survey. These stars demonstrate long-term velocity 
stability over a wide range of spectral types. Gl 694 (M2.5, V=10.5) 
shows that our velocity precision is reduced around faint, late spec- 
tral type stars, but that it is sufficient to detect the planet around 
Gl 179 (M3.5, V=11.96). The binned velocities with measurement 
uncertainties (but not jitter) are plotted. Panels are labeled with 
star name, spectral type, and rms to a linear fit. 

have median single measurement uncertainties of 1.2- 
1.4 ms~^. In some instances we made two or more con- 
secutive observations of the same star and binned the 
velocities in 2 hr time bins, thereby reducing the associ- 
ated measurement uncertainties. Gl 179 is substantially 
fainter = 11.96), and its median measurement uncer- 
tainty of 4.2 ms~^ is dominated by Poisson noise from 
photon statistics. For each of these stars, the measure- 
ment uncertainty is the weighted standard deviation of 
the mean of the velocity measured from each of the ^700 
2 A chunks in the echelle spectrum (Butler et al. 1996). 

4.2. HET Observations and Doppler Reduction 

We observed two of the stars, HD 34445 and G1179, 
with the Hobby-Ebcrly Telescope (HET) (Ramsey et al. 
1998) at McDonald Observatory, as part of the on-going 
Doppler surveys of FGK-type stars (Cochran et al. 2004) 
and M dwarfs (Endl et al. 2003). We use the High- 
Resolution-Spectrograph (HRS) (Tull 1998) in combina- 
tion with an 12-ccll at a resolving power of i? = 60, 000 
to measure precise radial velocities. The observing strat- 
egy and data reduction pipeline are detailed in Cochran 
et al. (2004) and a description of our l2-data modeling 



code "Austral" can be found in Endl et al. (2000). 

HET observations of HD 34445 begun in 2003 October 
and we collected a total of 45 HRS spectra since then. 
We started RV monitoring of Gl 179 a year later, in 2004 
October, and so far we have obtained 11 spectra of this 
faint M dwarf. 

4.3. Keplerian Models 

For each of the four stars below, we present a single- 
planet Keplerian fit to the measured velocities. For 
HD 126614 A, we add a linear velocity trend to the model. 
Table 3 summarizes the observing statistics, best-fit Kep- 
lerian orbital parameters, and measures of the goodness- 
of-fit for these planets. In this subsection, we describe 
the Keplerian modeling and related topics associated for 
all of the stars. We discuss the planet-bearing stars in- 
dividually in § 5-8. 

With the measured velocities for each star we per- 
formed a thorough search of the orbital parameter space 
for the best-fit, single-planet Keplerian orbital model us- 
ing the partially-linearized, least-squares fitting proce- 
dure described in Wright & Howard (2009). Each veloc- 
ity measurement was assigned a weight, w, constructed 
from the quadrature sum of the measurement uncertainty 
(aRv) and a jitter term (ajitter), i-e. w = 1/ (o-Rv+cr/itter)- 
Following Wright (2005), and based on the values of 
Shk, Mv, and B — V, we estimate ajitter for each of 
the stars (Table 3). These empirical estimates account 
for RV variability due to myriad sources of stellar noise, 
including the velocity fields of magnetically-controlled 
photospheric turbulence (granulation, super-granulation, 
and meso-granulation), rotational modulation of stellar 
surface features, stellar pulsation, and magnetic cycles, 
as well as undetected planets and uncorrected system- 
atic errors in the velocity reduction (Saar ct al. 1998; 
Wright 2005). We adopt larger values of jitter for pre- 
upgrade Keck measurements (crjitter = 3ms~^) than for 
HET measurements and post-upgrade Keck measure- 
ments (cTjittor = 2-2.5ms~^). The difference accounts for 
the slightly larger systematic errors incurred in the re- 
duction of prc-upgradc spectra. 

The Keplerian parameter uncertainties for each planet 
were derived using a Monte Carlo method (Marcy et al. 
2005) and do not account for correlations between pa- 
rameter errors. 

For stars with small Doppler amplitudes (HD 34445 
and HD 126614 A), we also explicitly considered the null 
hypothesis — that the velocity periodicity represented by 
the Keplerian orbital fit arose from chance fiuctuations 
in the velocities — by calculating a false alarm probabil- 
ity (FAP) based on Ax^, the goodness of fit statistic 
(Howard et al. 2009; Marcy et al. 2005; Gumming et al. 
2008). These FAPs compare the measured data to 1000 
scrambled data sets drawn randomly with replacement 
from the original measurements. (The linear velocity 
trends seen in HD 126614 A was subtracted from the ve- 
locities before scrambling.) For each data set we compare 
a best-fit Keplerian model to the null hypothesis (a linear 
fit to the data) by computing Ax^ =Xun "Xxep' where 
X^in and Xkcp values of x^ for linear and Ke- 

plerian fits to the data, respectively. The Ax^ statistic 
measures the improvement in the fit of a Keplerian model 
compared to a linear model of the same data. The FAP 
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TABLE 3 

Single Planet Keplerian Orbital Solutions 



r a,rciinctcr 








m 1 7Q 


P (yr) 


2.87 ±0.03 


3.41 ±0.05 


11.5±1.1 


6.26 ±0.16 


P(d) 


1049 ±11 


1244 ±17 


4218 ±388 


2288 ± 59 


e 


0.27 ±0.07 


0.41 ±0.10 


0.02 ±0.05 


0.21 ±0.08 


K (ms-l) 


15.7±1.4 


7.3 ±0.7 


23.3 ±1.4 


25.8 ±2.2 


Tp (JD - 2,440,000) 


13,781 ±48 


13,808 ±52 


14,494 ± 904 


15,140 ±104 


(deg) 


104 ± 19 


243 ± 19 


290 ± 78 


153 ±24 


dv/dt (m s"-*^ yr"-*^) 


= 0.0 


16.2±0.2 


= 0.0 


= 0.0 


Sim (^ivijupj 




n -1- n c\A 
u.oo m U-U-^: 




u-oz m u.u ( 


a (AU) 


2.07 ±0.02 


2.35 ±0.02 


5.15±0.29 


2.41 ±0.04 


A^obs (Keck, binned) 


68 


70 


39 


30 


Nohs (HET, binned) 


50 






14 


Median binned uncertainty (Keck, ms~^) 


1.4 


1.3 


1.3 


4.7 


Median binned uncertainty (HET, ms~^) 


2.7 






7.2 


Assumed Keck pre-upgrade jitter (ms~^) 


3.0 


3.0 


3.0 


3.0 


Assumed Keck post-upgrade jitter (ms~^) 


2.0 


2.0 


2.0 


2.5 


Assumed HET jitter (ms"-"^) 


2.0 






2.5 


rms to fit (ms^-*-) 


7.31 


3.99 


3.31 


9.51 




2.39 


1.42 


1.01 


1.78 



is the fraction of scrambled data sets that have a larger 
value of Ax^ than for the unscrambled data set. That 
is, the FAP measures the fraction of scrambled data sets 
where the statistical improvement from a best-fit Keple- 
rian model over a linear model is greater than the statisti- 
cal improvement of a Keplerian model over a linear model 
for the actual measured velocities. We use Ax^ as the 
goodness-of-fit statistic, instead of other measures such 
as Xv for ^ Keplerian fit, to account for the fact that the 
scrambled data sets, drawn from the original velocities 
with replacement, have different variances, which some- 
times artificially improve the fit quality (i.e. some scram- 
bled data sets contain fewer outlier velocities and have 
lower rms). It is important to note that this FAP does 
not measure the probability of non-planetary sources of 
true velocity variation masquerading as a planetary sig- 
nature. 

5. HD 34445 

5.1. Stellar Characteristics 

HD 34445 (= HIP 24681) is spectral type GO, with 
y = 7.31, B-y = 0.66, and Afy=4.04, placing it -0.8 
mag above the main sequence as defined by Wright 
(2005). It is chromospherically quiet with 5'hk = 0.15 
and log = —5.07 (Isaacson 2009), implying a rota- 
tion period of ~22 days (Noyes et al. 1984). Valenti 
& Fischer (2005) measured a super-solar metallicity of 
[Fe/H] = +0.14 using SME. Its placement above the main 
sequence and low chromospheric activity are consistent 
with an old, quiescent star of age 8.5 ± 2 Gyr. 

5.2. Photometry from Fairborn Observatory 

We have 560 individual photometric observations of 
HD 34445 spanning seven consecutive observing seasons. 
The mean short-term standard deviation of HD 34445 is 
comparable to the mean standard deviations of the com- 
parison stars (Table 2; columns 7 & 8); both are within 
the range of our typical measurement precision for sin- 
gle observations. We performed periodogram analyses 
over the range of 1-100 days on the seven individual ob- 
serving seasons and on the entire data sets and found no 



HD 34445 




RMSk.,1, = 6.40 ms'' 



1998 2000 2002 2004 2006 2008 2010 
Time (Years) 

Fig. 3. — Measured velocity vs. time for HD 34445 with the as- 
sociated best-fit Keplerian model {dashed line). Filled blue circles 
represent measurements from Keck, while filled red triangles rep- 
resent HET measurements. The error bars show the quadrature 
sum of measurement uncertainties and jitter. 

significant periodicity in either HD 34445 or the compari- 
son stars. Therefore, our short-term variability measure- 
ments of 0.00155 and 0.00156 mag are upper limits to 
any real nightly brightness variability in the target and 
comparison stars. 

Similarly, the long-term standard deviations of 
HD 34445 and the comparison stars are similar (Table 2, 
columns 9 & 10), indicating that we have not resolved 
any intrinsic year-to-year variability in HD 34445. Henry 
(1999) and Wright et al. (2008) present data sets that 
have seasonal means with standard deviations as low 
as 0.0002 mag, which we take as our mesurement pre- 
cision for seasonal mean magnitudes. Thus, our non- 
detection of long-term variability in HD 34445 may be 
compromised by low-level variability in one or both of 
the comparison stars. 

5.3. Doppler Observations and Keplerian Fit 

We monitored the radial velocity of HD 34445 at Keck 
for the past 12 yr and at the HET for the past 6 yr. Table 
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4 gives the Julian Dates of the observations, the mea- 
sured (relative) radial veloeities, the assoeiated measure- 
ment uncertainties (excluding jitter), and the telescope 
used. The 118 velocities listed in Table 4 have an rms 
of 12.0ms~^and are plotted as a time series in Figure 3 
along with the best-fit, single-planet Keplerian model. 

The velocities plotted in Figure 3 reveal a ~2.8yr 
periodicity that is apparent by visual inspection. We 
searched periods near 2.8 yr, and a wide variety of other 
periods, to find a single-planet Keplerian model with 
best-fit parameters, P = 2.87 ± 0.03 yr, e = 0.27 ±0.07, 
and iir=15.7±1.4ms~"'^, implying a planet of minimum 
mass Af sin z = 0.79 Mjup orbiting with semimajor axis 
a = 2.07AU. The full set of orbital parameters is listed 
in Table 3.^^ This joint fit has an rms of 7.31 ms~^ with 
Xv = 2.39. The planetary signal was detected in the 
Keck and HET data sets individually, which have an rms 
of 6.39 ms~^ and 7.95 ms~^ about the joint fit, respec- 
tively. 

We also considered the null hypothesis — that the ob- 
served RV measurements are the chance arrangement of 
random velocities masquerading as a coherent signal — by 
calculating an FAF. As described in §4.3, we computed 
the improvement in Ax^ from a constant velocity model 
to a Keplerian model (without a trend) for 10'' scram- 
bled data sets. We found that no scrambled data set had 
a larger value for Ax^ than for the measured velocities, 
implying an FAP of less than ^O.OOl for this scenario. 

While the single-planet model appears secure, the rms 
of the velocity residuals (6-8 ms~^) is a factor of ~2 
higher than expected based on the measurement uncer- 
tainties and jitter. Two possible explanations for this ex- 
cess variability are underestimated jitter and additional 
planets. Jitter seems an unlikely explanation given this 
star's metallicity, color, and modest evolution. As a 
comparison, the 5 stars most similar to HD 34445 {B - 
V and My within 0.05 mag) with 10 or more Keck ob- 
servations have velocity rms in the range 2.6-4.9 ms~^. 
We searched a wide range of possible periods for a sec- 
ond planet and found several candidates, the strongest 
of which has Pc = 117 d, sin i = 52 , and an FAP of 
a few percent. A significant number of additional mea- 
surements are required to confirm this planet and rule 
out other periods. 

6. HD 126614 A 

6.1. Stellar Characteristics 

HD 126614 (= HIP 70623) is identified in the Henry 
Draper and Hipparcos catalogs as a single star. As de- 
scribed in § 6.5, we directly detected by adaptive optics 
(AO) a faint M dwarf companion separated by 489 mas 
from the bright primary. We refer to the bright pri- 
mary star as HD 126614 A and the faint companion as 
HD 126614 B. These stars are unresolved in the Doppler 
and photometric variability observations described be- 
low. The planet announced below orbits HD 126614 A 
and is named HD 126614 Ab. In addition, HD 126614 A 
is orbited by a second M dwarf, NLTT 37349, in a much 
wider orbit (Gould & Chaname 2004). This outer stellar 
companion is separated from HD 126614 A by 42" and 

Advanced mention of the existence of this planet for HD 34445 
was made in Fischer & Valenti (2005). 
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Fig. 4. — Measured velocity vs. time for HD 126614 A (filled cir- 
cles) and best-fit Keplerian model {dashed line). The strong linear 
trend of dv/dt = l().2ms^^ yv^^ is clearly seen. The error bars 
that represent the quadrature sum of measurement uncertainties 
and jitter are not visible on this scale. These data are plotted again 
with the linear trend removed in Figure 5. 

docs not contaminate the Doppler or photometric obser- 
vations. 

HD126614A is spectral type KO, with V = S.81, B- 
F = 0.81, and My = 4.64, placing it ~1.2 mag above 
the main sequence. It is chromospherically quiet with 
•S'hk = 0.15 and logi?^j^ = -5.44 (Isaacson 2009), imply- 
ing a rotation period of 99 d, which is off the scale of 
the calibration, but suggests that it is longer than 50 d 
(Noyes et al. 1984). Valenti & Fischer (2005) measured 
an extremely high metallicity of [Fe/H] = -1-0.56 using 
SME. This is the highest metallicity measured in the 1040 
stars in the SPOCS catalog (Vaicnti & Fischer 2005). 
The low chromospheric activity (S'hk = 0.15) is consis- 
tent with the age estimate of 7.2 ± 2.0 Gyr from the SME 
analysis (Valenti & Fischer 2005). 

6.2. Photometry from Fairborn Observatory 

We have only 113 photometric observations of 
HD 126614 A covering two observing seasons. Peri- 
odogram analyses over the range of 1-100 days found no 
significant periodicity in cither HD 126614 A or its com- 
parison stars. The short-term variability of 0.00158 mag 
in the target star is comparable to 0.00163 mag in the 
comparison stars and to our measurement precision. The 
measurement of long-term variability in HD 126614 A is 
compromised by low-amplitude variability in one or both 
comparison stars. Thus, 0.0016 mag serves as an up- 
per limit to both long- and short-term variability in 
HD126624A. 

6.3. Doppler Observations and Keplerian Fit 

We monitored the radial velocity of HD 126614 A at 
Keck for the past 11 yr. Table 5 gives the Julian Dates 
of the observations, the measured (relative) radial ve- 
locities, and the associated measurement uncertainties 
(excluding jitter). The 70 velocities listed in Table 5 and 
plotted as a time scries in Figure 4 display a strong linear 
trend of 16.2ms^^ yr^^ with an rms of 6.6ms~^ about 
that trend. For clarity, the same velocities are plotted in 
Figure 5 with the trend subtracted. 
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Fig. 5. — The same measured velocities (filled cirvles) and single 
planet Keplerian model (dashed line) for HD 126614 A as plotted 
in Fig. 4, except the linear trend was subtracted before plotting. 

After subtracting the linear trend, a Lonib-Scargle pe- 
riodogram of the HD 126614 A velocities reveals a strong 
periodicity near 3.4 yr. This periodicity closely matches 
the period of the best-fit, single-planet Keplerian model, 
which was found after a thorough search of the Keple- 
rian parameters. Our model of a single planet plus a lin- 
ear velocity trend has a distinctive staircase-like appear- 
ance (Figure 4) and the following best-fit parameters: 
P = 3.41 ± 0.05 yr, e = 0.41±0.10, ii' = 7.3±0.7ms-i, 
and dv/dt = 16.2±0.2ms~^ yr~^ , implying a planet of 
minimum mass M sin i = 0.38 Mjup orbiting with semi- 
major axis a = 2.35AU. The full set of orbital parame- 
ters is listed in Table 3. The time series velocities (with 
trend subtracted) are plotted with the Keplerian model 
in Figure 5. 

This model of a single planet plus linear trend has rms 
velocity residuals of 3.99 ms"'^ and Xv = 1-42, implying 
an adequate fit. We computed an FAP for the single- 
planet model using the Ax^ statistic, as described in 
§ 4.3. We found that no scrambled data set had a larger 
value for Ax^ than for the measured velocities, implying 
an FAP of less than ~0.001 for this scenario. 

6.4. Source of Linear RV Trend 

The significant linear trend of 16.2 ms~^ yr^^ is likely 
due to a long-period stellar or planetary companion. 
We tried fitting the observed velocities with a 2-planet 
model, using the Keplerian parameters of HD 126614 Ab 
and an outer planet with a wide variety of long peri- 
ods as initial guesses. We were unable to find any two- 
planet models that improved Xv with statistical signifi- 
cance. Put another way, with a time baseline of 10 yr we 
do not detect curvature in the residuals to the one-planet 
model. 

We also considered NLTT 37349 (2MASS J14264583- 
0510194) as the companion responsible for the ob- 
served acceleration. Gould & Chaname (2004) iden- 
tified this object as a common proper-motion part- 
ner to HD 126614 A with a sky-projected separation of 
42". The observed photometric properties of this ob- 
ject are My = 12.02, V = 16.19, and F - J = 
4.00. Presumably NLTT 37349 has the same paral- 
lax (13.8mas) and distance (72.4pc) as HD126614A. 



Taken together, these properties are consistent with an 
M dwarf, roughly M4. implying a mass of m ^ 0.2 Mq. 
To compute the approximate gravitational influence of 
NLTT 37349, we assume that the line of sight separa- 
tion between HD 126614 A and NLTT 37349 is compara- 
ble to their physical separation in the sky plane, rscp ~ 
3000 AU. The acceleration predicted by this model, fV ^ 
Gm/rggp = 0.004 m s~^ yr~^, is much too small to ac- 
count for the observed acceleration of 16.2ms~^yr~^. 
Thus, NLTT 37349 is not the source of the observed lin- 
ear velocity trend. 

Secular acceleration is a potential cause of a linear ve- 
locity trend, especially for stars with significant proper 
motion (such as HD 126614 A). While our standard ve- 
locity pipeline removed this effect from the velocities re- 
ported Table 5, it is reassuring to confirm that the cal- 
culated secular acceleration is inconsistent with the ob- 
served velocity trend. (The pipeline also removed the 
motion of Keck Observatory about the barycentcr of the 
Solar System.) As discussed in Kiirster et al. (2003) 
and Wright & Howard (2009), a star's proper motion 
will cause the the radial component of its space veloc- 
ity vector to change with position on the sky, result- 
ing a secular acceleration of = DfP (to first order), 
where Vr is the bulk radial velocity of the star, D is 
the star's distance, and ii is the total proper motion in 
radians per unit time. For the most extreme cases of 
nearby stars with significant proper motion, fV can be 
as high as a few ms~"'^yr~^. For HD126614A, we find 
Vr ~ 0.1ms~^yr~^, ruling out secular acceleration as 
the cause of the observed 16.2 ms~^ yr~^ trend. 



6.5. AO Observations 

Having ruled out the above reasons for the linear RV 

trend, we considered a stellar companion close enough 
to HD126614A to have been missed by prior imaging 
surveys. To search for such a companion, we obtained 
direct imaging of HD 126614 A in J, H, and K-short (Kg) 
bands on the Hale 200" Telescope at Palomar Observa- 
tory on 13 April 2009 using the facility adaptive optics 
imager PHARO (Hay ward et al. 2001). In each band, 
we co-added approximately 150 431-ms exposures. As 
expected, the AO correction was best in Kg-band, trans- 
lating into the smallest errors. A nearby faint companion 
~500 mas to the Northeast of HD 126614 A was identified 
in each band by visual inspection (Figure 6). We con- 
structed an empirical PSF in each band from images of 
a calibrator, and images were fit with a two-PSF model. 
Parameter errors were calculated using the curvature of 
the surface, and for consistency, the error in contrast 
ratio compared to the percent flux of the fit residuals. 
Combining the astromctric fits in all bands, we find a 
projected separation of 489.0 ± 1.9 mas at a position an- 
gle of 56.1 ± 0.3 dcg. 

We are unaware of a previous detection of this 
faint companion to HD 126614 A, which we name 
HD 126614 B. We summarize the astrometric and pho- 
tometric results of the AO observations and the inferred 
properties of HD 126614 B in Table 6. 

HD 126614 A and HD126614B arc unresolved in the 
spectroscopic and photometric observations described in 
previous sub-sections. We estimate the amount of V- 
band contamination from HD 126614 B in these obser- 
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Fig. 6. — Direct observations with adaptive optics of HD 126614 A and HD 126614 B taken with the PHARO imager on Hale 200" Telescope 
at Palomar Observatory. The three panels show Ks-band images of the target image (left), the best two-PSF fit (center), and fit residuals 
(right). Each image is oriented North-up and East-left. The fainter star, HD 126614 B, is clearly detected to the Northeast of the brighter 
target star, HD 126614 A. The vertical and horizontal axes of each image are labeled with detector pixels (the plate scale is 25.2 mas/pixel). 
The image grcyscalc has been chosen to highlight the companion. Airy rings, and diffraction spikes; the residuals (right panel) have been 
stretched by a factor of 150. The full width half-maximum of the AO-corrected PSF was ~100 mas in each band (with Strehl ratios of 30% 
in Kg-band and 10% in J-band). 



TABLE 6 

Summary of AO imaging observations and isochrone analysis of HD 126614 B. 





HD 126614 B 






G1661B 


Padova Model 


Filter 


Separation 


Pos. Angle 


Cent. Ratio 


MjHK 


MjHK 


HD 126614 B Mass 




(mas) 


(deg) 




(mag) 


(mag) 


(Mq) 


J 


498.7 ± 6.1 


57.2 ± 0.9 


45.3 ± 4.1 


7.42 


7.46 


0.336 ± 0.015 


H 


487.1 ± 3.3 


56.1 ± 0.6 


38.8 ± 1.8 


6.97 


6.97 


0.307 ± 0.007 


K, 


488.5 ± 2.4 


56.0 ± 0.3 


35.0 ± 1.1 


6.75 


6.72 


0.308 ± 0.005 


Combined 


489.0 ± 1.9 


56.1 ± 0.3 








0.324 ± 0.004 



Note. — The absolute JHK magnitudes of HD 126614 B are derived from the measured contrast ratios, 
2MASS photometry, and the Hipparcos parallax of HD 126614 A. The .I-band results may be complicated 
by the presence of a diffraction spot near the position of the secondary. The uncertainties listed are based 
solely on the uncertainties in the contrast ratios. Accounting for all uncertainties, we estimate the combined 
mass to be accurate to perhaps 10%. For comparison, we show the absolute ,1HK magnitudes of G1661 B, 
with Mv — 11.3. This star was first identified in Reid & Gizis (1997) by speckle interferometry; we use 
updated photometry from Reid & Hawley (2005). 



vations using three methods: model isochrones, obser- 
vations of a similar star, and an empirical calibration 
employing metallicity effects. 

First, we interpolated the log(age) = 9.1, Z = 0.070 
(corresponding very nearly to [Fe/H] = -1-0.56) Padova 
isochrones (Girardi et al. 2000) at these absolute NIR 
magnitudes. These interpolations yield mass estimates 
for HD126614B (Table 6), the average of which is 
0.324 ± 0.004 M0, making HD 126614 B an M dwarf. The 
weighted average V-band magnitude implied by these 
measurements is My = 10.65 ± 0.03. 

Second, we find that the M dwarf Gl 661 B has nearly 
identical NIR absolute magnitudes to HD 126614 B (Ta- 
ble 6; Reid & Hawley 2005). Nonetheless, its absolute V 
magnitude is in severe disagreement with the Z = 0.070 
isochrones (and in even worse agreement with the solar 
metallicity isochrones). This suggests that the Padova 
isochrones significantly overestimate the V-band flux of 
M dwarfs. We thus consider the isochrone V magnitudes 
to be unreliable. 

Finally, we have applied the M dwarf photometric 
calibration of Johnson & Apps (2009) to the absolute 
Kg magnitude derived from the AO photometry. This 
empirical calibration is based on G-M binaries with 
well-measured metallicities (from LTE spectral synthe- 



sis analysis of the primary) and Kg-band magnitudes of 
both components. Johnson & Apps find that a star with 
[Fe/H] = +0.56 and Mk, = 6.75 should have My = 12.2. 
This estimate, which we feel is the most robust of the 
three we have performed, implies V = 16.5 and therefore 
negligible contamination in our optical spectra (AV = 
7.8, <0.1% contamination). 

Knowledge of the mass of HD 126614 B from the 
isochrones combined with the observed RV trend of 
HD 126614 A allows us to constrain the true physical sep- 
aration between the stars. Let 6 be the angle between 
the line of sight to the primary and the line connecting 
the primary to the secondary, such that the secondary 
is behind the primary when 9 — 0. The observed radial 
velocity trend (the instantaneous acceleration along the 
line of sight) is related to 9, the mass of the secondary 
Mb, and the true physical separation tab by 

• . ... 

t^A = — 2 COSt^, (1) 

''ab 

where tab is related to the apparent angular separation 
Pab through 
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and tta, the Hipparcos parallax of HD 126614 A. Given 
the observed radial acceleration of 16.2ms^^yr^^ and 
the best separation and mass from Table 6, there are two 
solutions for the implied physical separation between the 
stars: tab = 40^J and 50^3 AU (where the uncertainty is 
dominated by the uncertainty in the parallax) , which we 
can crudely combine as 45^g AU with the caveat that the 
probability distribution function is highly non-Gaussian. 

Holman & Wiegcrt (1999) studied the dynamical sta- 
bility of planetary orbits in binary systems. In the terms 
of that work, we find the mass ratio of the binary sys- 
tem to be /i = 0.22, at the edge of Hill's Regime. Hol- 
man & Wiegert find that for a binary in a circular or- 
bit, planets in S-type orbits of semimajor axis at are 
generally stable when secondary stars orbit with semi- 
major axes flB > a;,/0.38. For our system, this yields 
ae > 6.2 AU, consistent with our AO astrometry. In- 
deed, the HD 126614 system unconditionally passes this 
weak test for stability for ee < 0.6; beyond this regime 
more detailed knowledge of the binary orbit is required 
for further analysis. 

Pfahl & Muterspaugh (2006) discuss the prevalence 
of planets in binary systems such as HD 126614. They 
note that the literature contains a few planet detections 
in systems separated by ~ 20 AU, the most similar to 
HD 126614 being HD 196885 (Fischer et al. 2009, sub- 
mitted). Thus, HD 126614 Ab appears to be in a stable 
orbit in a binary system that is not atypical, even among 
systems with detected planets. 

6.6. Planetary Interpretation 

The identification of HD 126614 as a stellar binary with 
a separation of ~0'.'5 (unresolved by Keck/HIRES) raises 
the question of the origin of the periodic RV signal de- 
scribed in §6.3. We investigated the possibility that the 
RV variation was the result of distortions in the spectral 
line profiles due to contamination of the HD 126614 A 
spectrum by the much fainter HD 126614 B spectrum 
modulated by an orbital companion of its own. We 
present multiple lines of reasoning to argue against this 
alternative explanation and in favor of the Jovian-mass 
planet orbiging HD 126614 A described in § 6.3. 

From the AO observations in JHK bands, we es- 
timated the V-band flux of HD 126614 B in several 
ways with the most reliable method giving a flux ra- 
tio (A/B) of 7.8 mag. Thus in the iodine region of 
the spectrum, the lines from HD 126614 B are fainter 
by a factor of ^10'^. To produce a Doppler signal with 
K = 7ms-^ when diluted by HD 126614 A, the the signal 
from HD 126614 B would have an amplitude of approx- 
imately _fir = 7 km s~^ (implying a companion orbiting 
HD 126614 B with mass 0.25 Mq, an M dwarf). Yet, the 
stellar lines of HD 126614A, with va'mi =1.52kms^^, 
are narrower than the amplitude of the hypothetical 
_ft' = 7kms~^ Doppler signal. Thus, the stellar line pro- 
files of HD126614A and HD126614B would not be 
blended, but separated from each other for most of the 
3.41 yr orbit. This inconsistency casts serious doubt on 
the blend scenario as an explanation for the observed RV 
variation. 

Nevertheless, we investigated the blend scenario fur- 
ther with a spectral line bisector span test (Torres et al. 
2005; Queloz et al. 2001). We chose a region of the 



spectrum (6795-6865 A) to the red of the iodine re- 
gion for increased sensitivity to the spectral features of 
HDl 26614 B (an M dwarf) and a lack of telluric con- 
tamination. We cross-correlated each post-upgrade Keck 
spectrum against the solar spectrum using the National 
Solar Observatory solar atlas (Kurucz et al. 1984). Prom 
this representation of the average spectral line profile we 
computed the bisectors of the cross-corrleation peak, and 
as a measure of the line asymmetry we calculated the "bi- 
sector span" as the velocity difference between points se- 
lected near the top and bottom cross-correlation peak. If 
the velocities were the result of a blend, we would expect 
the line bisectors to vary in phase with the 3.41 yr period 
with an amplitude similar to K = 7ms^^. Instead, we 
find no significant correlation between the bisector spans 
and the RVs after subtracting the 16.2 ms~^ yr~^ trend. 
The Pearson linear correlation coefficient between these 
two quantities, r = — 0.28, demonstrates the lack of cor- 
relation. We conclude that the velocity variations are 
real and that the star is orbited by a Jovian planet. 

7. HD 13931 

7.1. Stellar Characteristics 

HD 13931 (= HIP 10626) is spectral type GO, with 
V = 7M, B-V = OM, and My = 4.32, placing it ~0.4 
mag above the main sequence. It is chromospherically 
quiet with 5'hk = 0.16 and log = —4.99 (Isaacson 
2009), implying a rotation period of ^26 days (Noyes 
et al. 1984). Valenti & Fischer (2005) measured an 
approximately solar metallicity of [Fe/H] = -1-0.03 using 
SME. Its placement above the main sequence and low 
chromospheric activity are consistent with an old, quies- 
cent star of age 6.4-10.4 Gyr. 

7.2. Photometry from Fairborn Observatory 

We have 181 photometric observations of HD 13931 
over two consecutive observing seasons. Periodogram 
analyses found no significant periodicity over the range 
1-100 days. The short-term standard deviations of the 
target and comparison stars (Table 2) provide an upper 
limit of 0.0016 mag for night-to-night variability. The 
long-term standard deviations suggest intrinsic low-level 
variability in HD 13931, but additional observations are 
needed to confirm it. 

7.3. Doppler Observations and Keplerian Fit 

We monitored the radial velocity of HD 13931 at Keck 
for the past 12 yr. Table 7 gives the Julian Dates of 
the observations, the measured (relative) radial veloci- 
ties, and the associated measurement uncertainties (ex- 
cluding jitter). The 39 velocities listed in Table 7 have 
an rms of 15.1 ms^"'^. These velocities arc plotted as a 
time series in Figure 7 along with the best-fit Keplerian 
model. 

A long period signal that has completed approxi- 
mately one cycle is clearly seen in the velocities plot- 
ted in Figure 7. After trying a wide variety of trial 
periods, we found the best-fit, single-planet Keple- 
rian model with P = 11.5 ± 1.1 yr, e 0.02 ± 0.08, and 
i4r = 23.3±2.9ms~^, implying a planet of minimum 
mass M sin i = 1.88 Mjup orbiting with semimajor axis 
a = 5.15AU. Note that the orbit is consistent with circu- 
lar. The full set of orbital parameters is listed in Table 
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Fig. 7. — Measured velocity vs. time for HD 13191 {filled circles) 
with the associated best-fit Keplerian model {dashed line). The 
error bars show the quadrature sum of measurement uncertainties 
and jitter. 

3. This model, with Xi^ = 1-01 and 3.31 ms^^ rms veloc- 
ity residuals, is an excellent fit to the data. Under the 
assumption of a one-planet model, we observed three RV 
extrema with Keck that strongly constrain the orbital 
parameters. 

8. GL179 

8.1. Stellar Characteristics 

G1179 HIP 22627) is spectral type M3.5, with 
y = 11.96, B-V=l.59, and My = 11.5, placing it -0.3 
mag above the main sequence. It is chromospherically 
active with S'hk = 0.96 (Isaacson 2009). Valenti & Fis- 
cher (2005) did not calculate stellar parameters using 
SME for G1179 because of its late spectral type. Simi- 
larly, the Noyes et al. (1984) calibration of stellar rotation 
does not apply tor B - V > 1.4. Johnson & Apps (2009) 
find a high metallicity of [Fe/H] = 0.30 ± 0.10 for G1179 
based on its absolute K-band magnitude, Mk, and V — K 
color. 

8.2. Photometry from Fairborn Observatory 

We have not made photometric observations of Gl 179 
from the APTs because at y = 11.96 the star is too faint. 

8.3. Doppler Observations and Keplerian Fit 

We monitored the radial velocity of Gl 179 at Keck for 
the past 10 yr and at the HET for the past 5yr. Table 
8 gives the Julian Dates of the observations, the mea- 
sured (relative) radial velocities, the associated measure- 
ment uncertainties (excluding jitter), and the telescope 
used. The 44 velocities listed in Table 8 have an rms of 
19.7 ms~^. Because G1179 is a faint M dwarf, the mea- 
surements have lower signal-to-noise ratios and larger un- 
certainties compared to the other stars in this paper. The 
measured velocities are plotted as a time series in Figure 
8 along with the best-fit Keplerian model. 

After trying a wide variety of trial orbital pe- 
riods, we find the best-fit, single-planet Keple- 
rian model has P = 6.26 ± 0.16 yr, e = 0.21 ±0.08, and 
X = 25.8 ± 2.2 ms~^, implying a planet of minimum 
mass M sin z = 0.82 Mjup orbiting with semimajor axis 
a = 2.41 AU. The full set of orbital parameters is listed in 
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Fig. 8. — Measured velocity vs. time for G1179 {filled circles) 
with the associated best-fit Keplerian model {dashed line). Filled 
blue circles represent measurements from Keck, while filled red 
triangles represent HET measurements. The error bars show the 
quadrature sum of measurement uncertainties and jitter. 

Table 3. We allowed for floating offsets between the HET 
RVs, the pre- upgrade Keck RVs, and the post-upgrade 
RVs when fitting the data. 

Note that the planet is clearly detected in the Keck RVs 
alone, but the HET measurements are crucial for the pre- 
cise determination of orbital parameters, especially ec- 
centricity and minimum mass. This single planet model, 
with Xi^ = l-78 and 9.51 ms^^ rms velocity residuals, is 
a statistically significant improvement over a model with 
no planets. However, the relatively large value of im- 
plies that either additional planets remain undetected in 
the system or we have underestimated the measurement 
uncertainties or jitter. We will continue to observe Gl 179 
to search for additional planets. 

9. DISCUSSION 

We present the detection of four new extrasolar plan- 
ets. These planets add to the statistical properties of 
giant planets in long-period orbits near the ice line. 
Long observational time baselines (10-12 yr) were nec- 
essary to accurately measure the high-eccentricity, low- 
amplitude signals of HD 34445b and HD 126614 Ab, as 
well as the intermediate-amplitude, but long-period sig- 
nals of HD 13931 b and Gl 179 b. 

HD 34445 b is a massive planet (M sin i = 0.79 Mjup) 
in a mildly eccentric (e = 0.27), long-period {P = 2.87 yr) 
orbit around an old GO star. We clearly detect this planet 
in the Keck and HET data sets individually, and their 
combination allows for more precise orbit determination. 
The relatively large residuals to the one-planet fit (rms — 
6-8 ms~^) hint at a second, unresolved planet in the sys- 
tem. Underestimated jitter would also explain the large 
residuals, but we deem this explanation less likely given 
the metallicity, color, and modest evolution of HD 34445. 

HD 126614 Ab is a massive planet {M sin z = 0.38 Mjup) 
in a long-period (P = 3.41yr), eccentric (e = 0.41) 
orbit around an extremely metal-rich star. At 
[Fe/H] = +0.56 ±0.04, HD 126614 A has the highest 
metallicity of the 1040 stars in the SPOCS catalog 
(Valenti & Fischer 2005). It also has the highest [Fe/H] 
of the ~250 stars with known planets and measured 
metallicities. We confirmed the high metallicity of 
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HD 126614 A by running a separate iodinc-frcc HIRES 
spectrum through the same SME pipchne used for the 
SPOCS catalog. We found [Fe/H] = +0.51 ±0.04, con- 
sistent with the SPOCS catalog value. Other authors 
also find an extremely high metallicity, including Castro 
et al. (1997) and Cenarro et al. (2007), who both find 
[Fc/H] = +0.55. These measurements, along with the de- 
tection of HD126614Ab, add statistical weight to the 
strong positive correlation between giant planet occur- 
rence and metallicity (Fischer & Valenti 2005). Indeed, 
HD 126614 A has been part of the planet-metallicity cor- 
relation story for some time. In an early paper dis- 
cussing the host star properties of some of the first cx- 
trasolar planets, Gonzalez et al. (1999) suggested that 
two bright, high-metallicity stars, namely HD 99109 and 
HD 126614, should be searched for Doppler variations. 
HD 99109 is known to host a planet with minimum mass 
Msini = O.SMjup and orbital period P=1.2yr (Butler 
ct al. 2006), and HD 126614 A now joins the list of high- 
metallity stars with planets. 

In addition to the planet orbiting HD 126614 A, we 
detected a faint M dwarf companion using adaptive 
optics and the PHARO camera at Palomar Observa- 
tory. This previously undiscovered star, HD 126614 B, 
has an estimated mass of 0.32 Mq and is separated 
from HD 126614 A by 489 ±1.9 mas at position angle 
56.1 ± 0.3 deg. This corresponds to a projected separa- 
tion of 33 AU. 

HD 13931 b is reminiscent of Jupiter in orbital pe- 
riod (P = 11.5yr), eccentricity (e = 0.02), and to a 
lesser extent mass (Msinz = 1.88Mjup). The host 
star, HD 13931, is also similar to the Sun in 
mass (M* = 1.02Mo) and metallicity ([Fe/H] = +0.03). 
HD 13931 b is one of only 4 known RV-detected plan- 
ets with orbital periods longer than 10 yr. The other 
such planets— 55 Cncd (Fischer et al. 2008), HD 217107c 
(Vogt et al. 2005; Wright et al. 2009) and HD 187123 c 
(Wright et al. 2007, 2009)— are all in multi-planet sys- 
tems. 

G1179b is a Jovian-mass (M sin i = 0.82 Mjup) planet 
in a long-period (P = 6.3yr) orbit. The host star, Gl 179, 
is one of only ^--^10 M dwarfs currently known to host a 
planet and is among the faintest (y = 11.96) stars with 
a planet discovered by RV measurements. This planet 
is detected in the Keck velocities alone, but without the 
HET measurements, the orbital parameters, especially 
eccentricity and minimum mass, would be determined 
more poorly. 

We note with interest that Gl 179 is an almost identical 
twin to G1876, an M dwarf known to host two Jovian 
planets locked in resonant orbits and a super-Earth in a 
P = 1.9d orbit (Marcy et al. 2001; Rivera et al. 2005). 
The stars are similar in effective temperature, mass, and 
age, as traced respectively by V—K color (5.00 for Gl 179 
and 5.15 for G1876), Mk (6.49 for G1179 and 6.67 for 
G1876_), and S'hk (0.96 for G1179 and 1.02 for G1876). 
The high metallicity of G1179 is also strikingly similar 
to the metallicity of G1876. Johnson & Apps (2009) 
estimate [Fe/H] = +0.3 and +0.37 for Gl 179 and G1876, 



respectively. Based on their analysis of M dwarfs with 
planets, Johnson & Apps find that the planet-metallicity 
correlation holds for stars at the bottom of the main- 
sequence. Gl 179 and its planet add statistical weight to 
this finding. 

The planet-bearing stars presented here are good can- 
didates for follow-up astrometric and direct imaging ob- 
servations. The astrometric perturbations from these 
planets on their host stars (lO's to lOO's of /ias) will 
be quite detectable with the expected sub-/ias sensitiv- 
ity of NASA's Space Interferometry Mission (SIM) or 
its proposed variants. Direct detection is also plausi- 
ble using coronagraphs/interferometers on space-borne 
and next-generation ground-based observatories, includ- 
ing GPI and SPHERE. HD 13931b may be the best can- 
didate with a maximum projected angular separation of 
120 mas. Indeed, SIM observations of HD 13931, even 
over a fraction of the orbital period (SIM has a planned 
mission duration of 5 years), combined with the RV mea- 
surements presented here, would completely determine 
the 3-dimensional orbit, giving accurate predictions of 
the best times for direct imaging observations. 
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TABLE 4 

Radial Velocities of HD 34445 



JD - 2440000 


R-adia-l Velocity 


Uncertctinty 


Telescope 




Cm s-i'l 


Cm 'i-l'l 




10838.76212 


-20.59 


1.83 


K 


11051.10569 


-1.84 


1.40 


K 


11073.04246 


-2.97 


1.21 


K 


11171.84795 


3.84 


1.55 


K 


11228.80594 


2.02 


1.47 


K 


11550.88478 


5.75 


1.35 


K 


11551.88096 


3.79 


1.43 


K 


11581.87017 


9.09 


1.54 


K 


11898.03487 


-20.63 


1.39 


K 


11974.76461 


-2.35 


1.53 


K 


12003.74221 


-13.51 


1.50 


K 


12007.72006 


-8.73 


1.41 


K 


12188.14432 


5.87 


1.59 


K 


12219.15109 


13.75 


1.79 


K 


12220.08113 


0.52 


1.80 


K 


12235.86269 


-3.92 


1.38 


K 


12238.88934 


3.14 


1.61 


K 


12242.92553 


6.49 


1.39 


K 


12572.99576 


21.87 


1.50 


K 


12651.93917 


2.68 


1.53 


K 


12899.09825 


-16.22 


1.29 


K 


12926.86832 


-23.82 


4.99 


H 


12932.00108 


-22.46 


2.09 


H 


12940.99257 


-18.74 


3.23 


H 


12942.98461 


-26.21 


2.37 


H 


12978.73611 


-24.64 


4.03 


H 


12979.73953 


-22.54 


4.01 


H 


12983.88276 


-13.97 


2.21 


H 


12984.88936 


-20.95 


3.41 


H 


12986.73101 


-6.99 


4.35 


H 


12986.86187 


-14.40 


2.99 


H 


13016.88898 


1.63 


1.21 


K 


13044.79245 


-6.83 


1.78 


K 


13255.96162 


14.40 


3.08 


H 


13258.95666 


12.95 


3.01 


H 


13260.96242 


15.60 


2.06 


H 


13262.94700 


12.35 


1.76 


H 


13286.88668 


8.19 


3.05 


H 


13303.12110 


11.45 


1.35 


K 


13338.90162 


-1.81 


1.21 


K 


13340.02404 


-3.95 


1.26 


K 


13359.69356 


8.66 


3.28 


H 


13365.68179 


12.38 


2.36 


H 


13368.95279 


5.50 


1.16 


K 


13369.80021 


2.74 


1.19 


K 


13756.75932 


7.60 


3.08 


H 


13775.69673 


-2.11 


4.85 


H 


13780.69018 


6.51 


2.08 


H 


13841.76181 


-0.41 


1.05 


K 


13978.97450 


-13.27 


2.42 


H 


13979.97609 


-19.30 


1.56 


H 


13982.10169 


-14.37 


1.16 


K 


13983.09819 


-8.77 


0.73 


K 


13984.10344 


-8.36 


0.89 


K 


1.3985.04721 


-10.67 


1.11 


K 


13985.98820 


-7.28 


3.29 


H 


13988.96142 


-12.49 


3.31 


H 


14003.92771 


-14.39 


1.91 


H 


14043.81148 


-23.50 


3.20 


H 


14055.77338 


-15.55 


2.73 


H 


14071.89776 


-6.16 


3.70 


H 


14096.67373 


-1.56 


3.40 


H 


14121.60247 


-0.97 


3.34 


H 


14130.87196 


-0.18 


1.17 


K 


14181.59907 


0.66 


5.45 


H 


14344.11108 


14.59 


1.16 


K 


14346.97829 


4.11 


1.86 


H 


14347.96596 


11.32 


3.26 


H 


14377.89952 


12.25 


1.64 


H 


14398.06266 


4.62 


1.17 


K 


14419.95414 


-5.91 


3.91 


H 


14452.85410 


-1.18 


2.39 


H 


14475.63706 


1.73 


3.76 


H 


14492. 8.'i216 


6,39 


1.25 


K 


14500.707:-!;-! 


3.81 


2.80 


11 
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TABLE 4 — Continued 



JD - 2440000 


Rciclitil Velocity 


Uncertainty 


TcIg S CO p 6 




(ms -*-) 


(Ulh j 




14544.81179 


17.87 


1.41 


K 


14546.73807 


19.12 


1.16 


K 


14547.60716 


19.93 


3.26 


H 


14730.92837 


4.20 


4.83 


H 


14762.00435 


3.81 


3.93 


H 


14777.95271 


4.68 


1.02 


K 


14781.79952 


7.11 


1.62 


H 


14791.00635 


15.34 


1.26 


K 


14807.92630 


-4.41 


1.46 


K 


14810.89575 


-0.73 


1.37 


K 


14838.84447 


-15.41 


1.40 


K 


14846.97512 


-11.01 


1.61 


K 


14856.59752 


-15.94 


2.22 


H 


14864.89937 


-19.53 


1.24 


K 


14865.83319 


-18.59 


0.96 


K 


14867.78335 


-21.11 


1.35 


K 


14877.68465 


-16.60 


2.26 


H 


14927.76208 


-12.33 


1.36 


K 


14928.72560 


-15.91 


1.23 


K 


14929.73178 


-7.89 


0.92 


K 


14934.73601 


-7.39 


1.36 


K 


15045.13099 


-11.01 


1.59 


K 


15049.12654 


-21.37 


1.43 


K 


15074.98600 


-24.11 


4.74 


H 


15076.12099 


-23.84 


2.69 


K 


15080.13683 


-25.10 


2.79 


K 


15084.13106 


-13.69 


2.79 


K 


15101.90535 


-34.62 


5.14 


H 


15110.12260 


-21.65 


1.46 


K 


15123.00538 


-6.66 


5.95 


H 


15135.10186 


-11.30 


1.37 


K 


15135.95876 


-10.74 


1.39 


K 


15142.80591 


-16.77 


4.28 


H 


15170.83820 


-8.05 


1.29 


K 


15187.83253 


-3.26 


1.27 


K 


1 Kl QQ QQ-1 CQ 

loioy.ooioy 


-D.UO 


l.O 1 


TT 
XV 


15190.67571 


-5.15 


4.38 


H 


15197.82141 


-11.12 


1.38 


K 


15215.60340 


-11.95 


2.12 


H 




-O.DO 




xV 


15231.95374 


-7.91 


1.21 


K 


15251.85907 


-6.17 


1.43 


is. 


15255.79845 


-3.31 


1.39 


K 




TABLE 5 






Radial Velocities of HD 126614 A 




JD - 2440000 


Radial Velocity 


Uncertainty 


Telescope 




(ms-l ) 


(ms-l ) 




11200.13355 


-175.78 


1.26 


K 


11311.92294 


-169.56 


2.17 


K 


11342.85107 


-155.44 


1.47 


K 


11370.81727 


-151.70 


1.44 


K 


11373.83647 


-146.70 


1.77 


K 


11552.15999 


-135.05 


1.28 


K 


11553.16885 


-136.56 


1.34 


K 


11581.17426 


-134.75 


1.26 


K 


11583.06953 


-140.19 


1.44 


K 


11585.12509 


-138.05 


1.07 


K 


11586.03518 


-134.02 


1.14 


K 


11680.01873 


-141.32 


1.54 


K 


11982.14743 


-126.57 


1.45 


K 


12003.94759 


-117.29 


1.48 


K 


12005.13470 


-118.83 


1.38 


K 


12065.93845 


-123.18 


1.49 


K 


12096.76464 


-116.69 


1.78 


K 


12128.75811 


-119.29 


1.23 


K 


12335.07744 


-112.03 


1.39 


K 


12363.08412 


-108.95 


1.52 


K 


12389.99741 


-110.09 


1.77 


K 


12445.84007 


-107.11 


1.37 


K 
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TABLE 5 — Continued 



JD - 2440000 


R-cidia.! Velocity 


Uncertaiinty 


Telescope 




(ms ) 


(Ulh ) 




12683.05902 


-82.87 


1.34 


K 


12711.99652 


-81.84 


1.40 


K 


12776.95616 


-76.98 


1.47 


K 


12849.80405 


-78.59 


1.51 


K 


13015.15707 


-70.77 


1.29 


K 


13016.16478 


-69.55 


1.20 


K 


13017.15230 


-70.34 


1.25 


K 


13179.83922 


-66.55 


1.47 


K 


13425.13398 


-57.71 


1.18 


K 


13777.11170 


-41.41 


1.02 


K 


13838.01968 


-39.62 


1.62 


K 


14131.15530 


-15.75 


1.24 


K 


14246.94875 


-18.96 


0.93 


K 


14247.94129 


-18.73 


1.21 


K 


14248.90277 


-18.55 


1.16 


K 


14251.83500 


-23.10 


1.13 


K 


14255.78231 


-20.33 


0.87 


K 


14277.75659 


-19.17 


1.13 


K 


14278.76238 


-18.63 


1.22 


K 


14285.77622 


-9.74 


1.28 


K 


14294.83121 


-18.26 


1.19 


K 


14300.80378 


-16.35 


1.15 


K 


14493.17109 


-9.46 


1.12 


K 


14549.01403 


-10.28 


1.26 


K 


14639.85723 


-11.05 


1.25 


K 


14839.12869 


-1.44 


0.87 


K 


14865.17003 


-3.63 


1.66 


K 


14866.08003 


-2.67 


0.91 


K 


14927.91558 


-3.71 


1.42 


K 


14929.09501 


-3.22 


1.31 


K 


14930.03891 


-2.54 


1.15 


K 


14964.02038 


-2.01 


0.73 


K 


14983.79390 


1.65 


1.23 


K 


14984.84861 


4.14 


1.19 


K 


14985.95010 


2.85 


1.38 


K 


14986.93908 


7.48 


1.23 


K 


14987.94064 


4.43 


1.30 


K 


14988.95513 


1.50 


1.28 


K 


15014.82834 


3.02 


1.24 


K 


15015.83799 


-1.02 


1.25 


K 


lOUiO-OO 1 Do 


n Qi 




TT 
IV 


15041.84683 


4.80 


1.33 


K 


15042.87063 


8.59 


1.45 


K 


15043.81939 


3.67 


1.33 


K 


iOU44. /yuio 




1.4U 


rv 


15190.15992 


19.99 


1.34 


K 


15197.15781 


21.27 


1.23 


TV- 
IV 


15229.04986 


27.36 


1.19 


K 




TABLE 7 






Radial Velocities of HD 13931 




JD - 2440000 


Radial Velocity 


Uncertainty 


'"Ppl p G f o n p 

J- v^l ^ 




(ms-l ) 


(ms-i ) 




10837.81679 


2.15 


1.59 


K 


11044.11226 


9.86 


1.26 


K 


11051.05061 


9.06 


1.27 


K 


11071.05334 


4.87 


1.66 


K 


11172.82617 


-1.57 


1.34 


K 


11374.13082 


8.45 


1.55 


K 


11412.08069 


2.92 


1.37 


K 


11543.80646 


1.85 


1.38 


K 


11552.79058 


-1.78 


1.56 


K 


11585.71983 


-5.01 


1.46 


K 


11899.88418 


-5.72 


1.33 


K 


12133.12439 


-13.19 


1.56 


K 


12236.81830 


-17.61 


1.44 


K 


12515.95280 


-25.82 


1.44 


K 


12574.90240 


-34.38 


1.45 


K 


12836.11950 


-35.49 


1.45 


K 


12989.72665 


-42.44 


1.35 


K 
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TABLE 7 — Continued 



JD - 2440000 


Rciditil Velocity 


Uncertainty 


TcIgS CO p 6 




I'm i;-! 


ycai, j 




13240.00935 


-36.93 


1.11 


K 


13961.11345 


-24.51 


0.64 


K 


14083.90911 


-24.85 


1.11 


K 


14397.92588 


-10.40 


1.05 


K 


14719.07441 


-6.08 


1.09 


K 


14809.83689 


-0.23 


1.14 


K 


14838.92165 


1.05 


1.08 


K 


14864.76650 


-1.77 


1.09 


K 


15015.10431 


5.07 


1.12 


K 


15016.10240 


3.72 


1.11 


K 


15017.10744 


3.23 


1.19 


K 


15042.12880 


8.88 


1.19 


K 


15044.13700 


8.69 


1.14 


K 


15049.04212 


9.74 


1.17 


K 


lOU 1 ( • iooUo 






XV 


15085.13348 


5.48 


1.02 


K 


15109.99945 


9.62 


1.38 


K 


15133.95126 


3.42 


1.25 


K 




fill 

D.ii 


L.ZZ 


IV 


15171.97729 


9.50 


1.59 


K 


i5io8. 89679 


8.46 


1.18 


TV- 
IV 


15231.78501 


6.06 


1.12 


K 




TABLE i 


i 






Radial Velocities of Gl179 




JD - 2440000 


Radial Velocity 


Uncertainty 






fms"-'- ) 


(ms-l ) 




11580.83131 


38.12 


4.05 


K 


11882.88790 


37.65 


4.87 


K 


11901.00250 


30.81 


5.13 


K 


12235.84867 


54.65 


6.14 


K 


12536.08795 


4.11 


4.50 


K 


12572.99093 


28.53 


5.14 


K 


12573.95028 


17.59 


4.52 


K 


12575.04686 


5.25 


4.48 


K 


12575.99081 


9.58 


5.42 


K 


12898.11579 


-8.75 


4.24 


K 


13014.81847 


-7.72 


4.53 


K 


13015.83165 


-16.60 


4.26 


K 


13016.83228 


-29.20 


3.72 


K 


13297.85507 


-1.09 


10.33 


H 


13299.84862 


14.31 


5.54 


H 


13302.97472 


2.28 


1.93 


K 


13314.80754 


8.21 


4.36 


H 


13340.72400 


2.54 


5.09 


H 


13631.92988 


10.84 


15.59 


H 


13644.90017 


25.11 


12.87 


H 


13686.92171 


9.61 


6.64 


H 


13719.83142 


28.29 


12.94 


H 


13984.08891 


32.57 


2.23 


K 


14130.85314 


35.38 


2.51 


K 


14397.93848 


39.52 


2.15 


K 


14411.93712 


17.96 


14.13 


H 


14771.80303 


22.44 


4.48 


H 


14778.99120 


21.36 


2.43 


K 


14790.99549 


13.51 


2.39 


K 


14807.91708 


12.78 


2.17 


K 


14838.99126 


19.51 


1.76 


K 


14846.95665 


11.56 


2.70 


K 


14864.95689 


14.55 


2.42 


K 


14928.73249 


6.67 


2.03 


K 


14929.72575 


4.21 


2.64 


K 


14934.73136 


-8.63 


2.88 


K 


15077.10989 


3.80 


1.99 


K 


15084.94943 


-1.33 


7.22 


H 


15170.78663 


-5.85 


1.68 


K 


15174.09252 


-14.84 


2.27 


K 


15175.70740 


-15.10 


8.34 


H 


15180.83537 


-32.31 


5.00 


H 


15187.83709 


-19.44 


2.46 


K 
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TABLE 8 — Continued 



JD - 2440000 


Radial Velocity 


Uncertainty 


Telescope 




(ms"-*- ) 


(ms-l ) 




15206.77569 


-3.89 


4.39 


H 



